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Abstract 

The recoil properties of the product nuclei from the interaction of 3.65 GeV/nucleon protons 
and deuterons from the Nuclotron and Synchrophasotron of the Laboratory of High Energies 
(LHE), Joint Institute for Nuclear Research (JINR) at Dubna with a ^^^Sn target have been 
studied using catcher foils. The experimental data were analyzed using the mathematical 
formalism of the standard two-step vector model. The results for protons are compared with 
those for deuterons. Our experimental results were compared to three different Los Alamos 
versions of the Quark-Gluon String Model (LAQGSM). The forward velocity v and the recoil 
nuclei kinetic energy increases linearly with increasing mass loss of the target AA, but seems 
to change its slope at around /S.A = 60. It seems that light- and medium-mass products are 
produced partly by a fragmentation mechanism. 

1 Introduction 

The investigation of interactions of high-energy particles with complex nuclei by the induced 
activity method is limited to measurements on the large residual nuclei that remain at the end 
of reactions. The study of large residual nuclei usually involves measurements of either their 
excitation functions or their recoil properties. 

In order to determine the recoil properties of nuclei "thick-target thick-catcher" experimens 
are used. In such experiments, the thicknesses of the target and catcher foils are larger than 
the longest recoil range. The quantities measured are the fractions F and B of product nuclei 
that recoil out of the target foil into the forward and backward directions, respectively. 

The results of the experiment are usually proceeded by the standard two-step vector repre- 
sentation The following assumptions are made in this model: 

(1) In the first step, the incident particle interacts with the target nucleus to form an 
excited nucleus with velocity f , momentum P, and excitation energy E* . 

(2) In the second step, the excited nucleus loses mass and excitation energy to form the 
final recoiling nucleus with an additional velocity V , which in general will have a distribution 
of values and directions. 

Usually, additional assumptions are made in most experiments: 

(a) The quantities v and P in the first step are constant 4] and lie in the forward 

direction. 
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(b) The velocity in the second step is isotropic. 
The results of the recoil experiments depend on the range-energy relation of the recoiling 
nuclei. It is convenient to express this relation as [2]: 



R = kV, (1) 

where R is the mean range (corresponding to V) of the recoil in the target material, k and n are 
constants and can be evaluated from tables of ranges of nuclei recoiling into various materials 
jH]. The following relations are used for the forward and backward fractions: 

FW = ^R[l + ^(n + 2)r/ + ^{n + 1) V]; BW = ^R[l - ^(n + 2)r] + ^{n + 1) V] (2) 

where r] = v/V and W is the target thickness in mg/cm?. 

The reaction product mean ranges [R) and the velocities transferred to residuals on the 
first (f) and second iV) steps of reaction are calculated using the following expressions 0: 

1 - |(n + 2)r7 + i(n + Ij^i]^ 4 

The recoil properties of non-fission reactions induced by protons with energy 1 GeV and 
above were analyzed with the two- vector model in Ref. Systematic deviations of the re- 
sulting parameters were attributed to the presence of fragmentation. The kinematic properties 
of radionuclides formed in photospallation reactions on complex nuclei at intermediate energy 
and a comparison with proton- nuclear reactions were made in Ref. |7j. We studied the recoil 
properties of nuclei produced in the photospallation of ^^Cu in Ref. f^. The purpose of the 
present work is to investigate the kinematic properties of product nuclei formed in the target 
^^^Sn bombarded with 3.65 GeV/nucleon protons and deuterons using catcher foils. 



2 Experimental setup and results 

Targets of enriched tin isotope ^^'^Sn were irradiated at the Nuclotron and Synchrophasotron 
of the LHE, JINR by proton and deuteron beams with energies of 3.65 GeV/nucleon. Irra- 
diations were of 6.42 hour for the proton beam and 1.083 hour for the deuteron beam. The 
deuteron beam had an elliptic form with axes of 3 and 2 cm. The proton beam had round form 
with a diameter of 2 cm. For beam monitoring, we employ the reactions ^''Al((i, 3p2n)^^Na 
and ^^Al(p, 3pn)^^Na whose cross sections are taken as of 14.2 ± 0.2mb |9j and 10.6 ± 0.8mb 
|10j , respectively. From the monitoring reactions the following beam intensities were obtained: 
0.768 X 10^^ d/hour and 0.114 x 10^^ p/hour. The total beam fluences were 3.21 x 10^^ protons 
and 2 x 10^^ deuterons. 

The target consisted of a high-purity target metal foil of size 20x20 mm^ sandwiched exactly 
by one pair of Mylar foils of the same size, which collected the recoil nuclei in the forward or 
backward directions with respect to the beam. The enrichment of the target was 98.7 %, the 
thickness of each target foil was 66.7 mg/cm? and the number of target piles was 15. The whole 
stack, together with an Al beam-monitor foil of 140 mg/ cm^ thickness was mounted on a target 
holder and irradiated in air. 

After irradiation the target foils and all of the forward and backward catcher foils from 
one target pile were collected separately, and assayed for radioactivities nondestructively with 



high-purity Ge detectors at LNP, JINR for one year. The radioactive nuclei were identified by 
characteristic 7 fines and by tfieir fialf-fives. Tfie spectra were evaluated witfi tfie code package 
DEIMOS32 [HI. 

Tfie kinematic cfiaracteristics of tfiirty product nuclei were obtained for deuteron- and 
proton-induced reactions. Tfie relative quantities of tfie forward- and backward-emitted nu- 
clei (relative to tfie beam direction) were calculated from relations: 

F = Np/iNt + Np + NB); B = NB/iNt + Np + Nb) (4) 

wfiere Np , Nb , Nt are tfie numbers of nuclei emitted in forward and backward catcfiers and 
formed in target foils, respectively. Tfie recoil parameters obtained in tfiese experiments are 
tfie forward-to-backward ratio, F/B, and tfie mean range, 2W{F + B). [Tfie mean range of 
tfie recoils is somewfiat smaller tfian 2W{F + B), but it is conventional to refer to tfie latter 
quantity as range]. Tfie matfiematical formalism of tfie standard two-step vector model j2] was 
used to proceed tfie experimental results. Tfie parameters k and n in equation (1) are obtained 
by fitting tfie range dependence on energy of accelerated ions witfiin tfie region from 0.025 to 5 
MeV/nucleon ^21- It is possible to calculate 1] and v from tfie equation (3), knowing tfie F/B 
ratio of tfie experiment. 

Our experimental results are sfiown in Tables 1 and 2. We note tfiat uncertainties concerning 
definite quantities in our tables are not listed to keep tfie tables concise. Tfiese uncertainties 
are about 10-15%. As is seen from tfie tables, tfie kinetic energies of product nuclei from 
proton-induced reactions are fiigfier tfian are tfie ones induced by deuterons. Probably, protons 
are more effective agents of linear momentum transfer on a per-nucleon basis wfien compared 
witfi deuterons [131 Ej. 

As sfiown in Fig. 1, tfie ratios F/B for botfi proton- and deuteron-induced reactions are of 
tfie order of ~ 3 — 4 for fieavy product nuclei and decrease to about ~ 2 for ligfit residuals. 
Sucfi a dependence could be explained by different mecfianisms for tfie production of nuclei in 
different mass regions. Ligfit nuclei may be produced by multifragmentation and evaporation 
tfiat lead to an isotropic distribution in tfie frame of excited residual nuclei, wfiile fieavy nuclei 
are produced mainly via tfie spallation mecfianism, witfi its products more in tfie forward 
direction. Some of our recent studies IT^ point to tfie multifragmentation mecfianism in 
tfie formation of product nuclei witfi ligfit and medium mass numbers. 

Our experimental results are compared witfi tfieoretical calculations by tfie LAQGSM03.01 
[HI, LAQGSM03.S1 [18j, and LAQGSM03.G1 (TH! models. 

LAQGSM03.01 [TTj is tfie latest modification of tfie Los Alamos version of tfie Quark-Gluon 
String Model JH] , wfiicfi in its turn is an improvement of tfie Quark-Gluon String Model |2D1 • It 
describes reactions induced by botfi particles and nuclei as a tfiree-stage process: IntraNuclear 
Cascade (INC), followed by preequilibrium emission of particles during tfie equilibration of 
tfie excited residual nuclei formed during tfie INC, followed by evaporation of particles from 
or fission of tfie compound nuclei. Tfie INC stage of reactions is described witfi a recently 
improved version ^7] of tfie time-depending intranuclear cascade model developed initially 
at Dubna, often referred in tfie literature simply as tfie Dubna intranuclear Cascade Model, 
DCM (see j2I] and references tfierein). Tfie preequilibrium part of reactions is described witfi an 
improved version of tfie Modified Exciton Model (MEM) from tfie Cascade-Exciton Model, 
GEM j2S|- Tfie evaporation and fission stages of reactions are calculated witfi an updated and 
improved version of tfie Generalized Evaporation Model code GEM2 by Furifiata wfiicfi 
considers evaporation of up to 66 types of different particles and ligfit fragments (up to ^^Mg). If 
tfie excited residual nucleus produced after tfie INC fias a mass number A < 11, LAQGSM03.01 



uses a recently updated and improved version of the Fermi Break-up model (in comparison with 
the version described in [201) calculate its decay instead of considering a preequilibrium stage 
followed by evaporation from compound nuclei, as described above. LAQGSM03.01 considers 
also coalescence of complex particles up to ^He from energetic nucleons emitted during the INC, 
using an updated coalescence model in comparison with the version described in 

LAQGSM03.S1 [l8j is exactly the same as LAQGSM03.01, but considers also multifrag- 
mentation of excited nuclei produced after the preequilibrium stage of reactions, when their 
excitation energy is above 2A MeV, using the Statistical Multifragmentation Model (SMM) by 
Botvina et al. fh] (the "S" in the extension of LAQGSM03.S1 stands for SMM). 

LAQGSM03.G1 ^18j is exactly the same as LAQGSM03.01, but uses the fission-like binary- 
decay model GEMINI of Charity et al. [211 5 which considers evaporation of all possible frag- 
ments, instead of using the GEM2 model PI (the "G" stands for GEMINI). 

As can be seen from Fig. 1, there is some disagreement between experimental data and theo- 
retical results by all three versions of LAQGSM considered here. In making such a comparison, 
we first recognize that the experiment and the calculations differ in that: 1) the experimental 
data were extracted assuming the "two-step vector model" pQ-jni, while the LAQGSM calcula- 
tions were done without the assumptions of this model; 2) the measurements were performed on 
foils (thick targets), while the calculations were done for interactions of protons/deuterons with 
nuclei (thin targets). These differences must be considered before assessing possible deficiencies 
in the models. 

It is interesting to note that the measured forward velocity v increases practically linearly 
with the increase of AA {AA = Atarg — Ares, where Atarg is the mass number of target and 
mass number of product nuclei), but seems to change it slope at around AA = 60. 
Comparison of these results with theoretical calculations shows that this is due to different 
mechanisms for the production of the measured nuclei (Fig. 2). 

Fig. 3 shows the dependence of the fragment kinetic energy, T^in, on the fractional mass 
loss AA/A. Gumming and Bachmann |2ZI and Winsberg P] have shown that Tkm should 
increase linearly with AA/A for reactions in which the velocity of the product is due to the 
vectorial addition of the randomly directed recoil velocities resulting from particle emission. 
For comparison. Fig. 3 shows also data measured for a Ag target by Porile et al. (see Table II 
in [0] and references therein). One can see that the Ag results agree well with our current ^^^Sn 
data. 

The deviation from a linear trend seen in Fig. 3 for large fractional mass losses probably 
indicates a change in mechanism of the production of light nuclides. The formation of light 
fragments from highly-excited nuclei is of a permanent interest in the literature and usually 
considers a multi-body breakup j2Hl- One possible mechanism for such a process would be a 
simultaneous clustering of nucleons into fragments near the liquid-gas critical point [22] • This 
process is essentially different from the sequential evaporation process by which deep spallation 
products are formed. It appears that the mean kinetic energies of the products provide a 
qualitative method for distinguishing between these two mechanism [2^. 

As one can see from Fig. 3, the results by LAQGSM03.S1, which considers multifragmenta- 
tion [25^ of excited nuclei when their excitation energy is above 2 MeV/nucleon, overestimate 
significantly the values of the measured mean kinetic energies of light products. This could 
be an indication that multifragmentation become important only at higher excitation energies, 
4-5 MeV/nucleon instead of 2 MeV/nucleon as considered by LAQGSM03.S1, in complete 
agreement with the very recent ISiS measurements [3G\ . 

The evaluated "experimental" excitation energies of residual nuclei produced after the first. 



cascade stage of reactions are shown in the last column of Tables 1 and 2. The relation between 
the excitation energy (E*) and v may be estimated as 

E* = 3.253 * 10-^k'ATv[Tp/{Tp + 2)]°■^ (5) 

where E* and the bombarding energy Tp are expressed in terms of irtpC^. At is the target mass 
in amu and v is in units of (MeV/amM)°'^. Usually, the constant k' is taken as /c = 0.8. 

As is seen from tables, the excitation energies estimated according to Eq. (5) are higher than 
the multifragmentation threshold Eth = 2 - A MeV/nucleon {Etot = 216 - 424 MeV) [31^ for 
the light and medium products. This could be an indication that light and medium fragments 
are produced not only via the evaporation mechanism but also via multifragmentation. 
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Table 1. Kinematic characteristics of product nuclei from deuteron-induced reactions 



Product 


/ 


n 


2W{F + B) 


Ti.„, (MeV) 




E* (MeV) 




2.03 


0.175 


4.08±0.40 


15.75± 2.37 


0.2255 


619.62 


28Mff 


1.77 


0.142 


4.49±0.70 


19.70± 4.80 


0.1913 


522.36 


42 


2.36 


0.212 


2.49±0.60 


12.07± 5.11 


0.1861 


496.32 




2.73 


0.246 


1.93±0.35 


7.54± 2.42 


0.1680 


449.30 




2.00 


0.172 


2.49±0.58 


12.21±4.63 


0.1457 


394.77 


529^11 


1.97 


0.168 


0.94±0.21 


2.86 ±1.03 


0.0678 


172.55 


56 Mn 


2.39 


0.215 


1.67±0.61 


6.91 ±4.06 


0.1212 


329.11 


'^^Ga 


2.18 


0.193 


1.51±0.44 


6.89 ±3.21 


0.0996 


270.13 


^^Se 


3.08 


0.274 


0.72±0.16 


2.25± 0.83 


0.0769 


210.08 




3.02 


0.269 


0.91±0.28 


3.16± 1.70 


0.0927 


252.67 




1.97 


0.167 


1.41±0.39 


6.51± 2.88 


0.0785 


212.61 




4.09 


0.338 


1.14±0.14 


4.85± 0.98 


0.1318 


361.96 


71 As 


2.84 


0.255 


1.45±0.27 


6.74± 2.03 


0.1261 


343.88 




3.39 


0.296 


1.11±0.21 


4.73± 1.45 


0.1130 


309.10 


85y 


3.51 


0.304 


0.89±0.12 


3.35± 0.72 


0.0963 


263.75 


86mY 


3.09 


0.275 


0.97±0.13 


3.88± 0.82 


0.0934 


254.93 


87mY 


3.44 


0.299 


0.97±0.17 


3.88± 1.12 


0.1010 


276.46 


86 Zr 


4.78 


0.372 


0.78±0.11 


2.78± 0.66 


0.1061 


292.67 


892r 


2.96 


0.265 


0.88±0.16 


3.33±0.97 


0.0821 


224.13 


90Mo 


2.39 


0.214 


0.84±0.17 


3.26±1.08 


0.0657 


178.49 


93m]yfQ 


2.07 


0.179 


0.52±0.07 


1.46± 0.33 


0.0362 


98.25 




3.28 


0.288 


0.45±0.09 


2.38± 0.48 


0.0515 


140.56 


s^Tc 


3.12 


0.277 


0.36±0.09 


0.86± 0.34 


0.0423 


115.53 


95Tc 


3.82 


0.323 


0.54±0.08 


1.61± 0.39 


0.0669 


183.76 


96Tc 


2.06 


0.179 


0.36±0.06 


0.85± 0.21 


0.0270 


73.29 


'^^Ru 


3.64 


0.312 


0.66±0.16 


2.28± 0.89 


0.0764 


209.48 


Rli 


3.56 


0.307 


0.15±0.07 


1.27± 0.33 


0.0555 


152.21 


104 


4.46 


0.357 


0.23±0.05 


0.45± 0.14 


0.0374 


103.03 


109jj^ 


3.24 


0.286 


0.26±0.04 


0.58± 0.14 


0.0332 


90.91 




2.00 


0.180 


0.21±0.05 


0.41± 0.16 


0.0168 


45.73 


iiij^ 


2.33 


0.210 


0.16±0.04 


0.26± 0.11 


0.0161 


43.80 




2.11 


0.180 


0.25±0.05 


0.56± 0.18 


0.0064 


56.20 



Table 2. Kinematic characteristics of product nuclei from proton-induced reactions 



Product 


/ 


n 


2W{F + B) 


Tkir, (MeV) 




E* (MeV) 




1.95 


0.165 


5.37± 0.35 


23.91± 2.36 


0.2628 


583.06 


28Mff 


2.03 


0.175 


4.61± 0.30 


20.54± 2.11 


0.2399 


530.01 




2.15 


0.189 


3.12± 0.20 


17.88± 2.04 


0.2023 


435.23 




2.19 


0.194 


2.37± 0.14 


10.83± 1.16 


0.1595 


343.19 




2.50 


0.225 


2.83± 0.18 


14.97± 1.57 


0.2111 


463.51 


46 Sc 


1.94 


0.165 


2.94± 0.20 


15.60± 1.68 


0.1547 


338.71 


^^Sc 


2.08 


0.181 


2.55± 0.34 


12.07± 2.60 


0.1467 


321.13 


48Y 


2.21 


0.196 


2.56± 0.19 


13.67± 1.67 


0.1687 


369.73 


^29Mn 


2.21 


0.196 


2.09± 0.14 


10.43± 1.09 


0.1415 


310.03 




3.21 


0.281 


2..38± 0.18 


12.66± 1.51 


0.2122 


167.95 


67Ga 


2.72 


0.245 


1.66± 0.11 


8.04± 0.84 


0.1362 


299.48 


^^As 


3.83 


0.324 


1.54± 0.10 


7.39± 0.79 


0.1665 


368.74 


7^Se 


3.96 


0.331 


0.97± 0.06 


3.62± 0.37 


0.1174 


260.25 


^^Se 


5.17 


0.389 


0.98=b 0.02 


3.77± 0.15 


0.1378 


307.88 


^^Br 
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Figure 1: F/B versus the fractional mass losses AA/At. a) for deuteron-induced reactions b) 
for proton-induced reactions. Solid lines (1) show calculations by LAQGSM03.01, dashed lines 
(2) by LAQGSM03.S1, and dotted lines (3) by LAQGSM03.G1. 
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Figure 2: Dependence of forward velocity on the number of emitted nucleons: a) for 
deuteron-induced reactions b) for proton-induced reactions. Solid lines (1) show calculations by 
LAQGSM03.01, dashed lines (2) by LAQGSM03.S1, and dotted lines (3) by LAQGSM03.G1. 




Figure 3: Dependence of the kinetic energy of the product nuclei on the fractional mass 
losses l^AjAt. a) for deuteron-induced reactions b) for the proton-induced reactions (•). For 
comparison, ■ show experimental results for the Ag target tabulated in Ref. 0. Solid lines (1) 
show calculations by LAQGSM03.01, dashed lines (2) by LAQGSM03.S1, and dotted lines (3) 
by LAQGSM03.G1. 



